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Screening masses of hot SU(N) gauge theory, defined as poles of the corresponding propagators are 
studied in 3d adjoint Higgs model, considered as an effective theory of QCD, using coupled gap 
equations and lattice Monte-Carlo simulations (for N = 2). In so-called A gauges non-perturbative 
evidence is given for the gauge independence of pole masses within this class of gauges. Application 
of screening masses in a novel resummation prescription of the free energy density is discussed. 
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I. INTODUCTION 

Finite temperature SU(N) theory undergoes a phase 
transition at some temperature T c from the confined to 
the deconfined phase. Above this temperature chromo- 
electric fields are screened with finite inverse screening 
length, the inverse of the so-called Debye mass. Al- 
though there is no confinement for temperatures above 
T c , naive perturbation theory is known to fail because 
of severe IR divergencies. The by-now well known re- 
summation techniques, though solve a part of this prob- 
lem (e.g. a weak coupling expansion can be obtained up 
g 5 order for the free energy) the resulting resummed se- 
ries shows very bad convergence properties |]|| (]. For 
the Debye mass IR problems imply that the naive defi- 
nition IIi(fco = 0, k — * 0) (where 11^ is the longitudinal 
self-energy ) is not applicable beyond the leading order. 
Rebhan has suggested to define the Debye mass as the 
pole of the logitudinal part of the propagator |Q. This 
definition yields gauge invariant results, however, it re- 
quires the introduction of the so-called magnetic mass, a 
concept, which was introduced long ago to cure IR diver- 
gencies due to static magnetic fields ||. Analogously to 
the electric mass the magnetic mass can be defined as the 
pole of the transverse gauge boson propagator. Although 
the magnetic mass is non-calculable in perturbation the- 
ory, numerical lattice studies of the finite temperature 
gluon propagator indicate its existence ||J^]. Also a self- 
consistent resummation of perturbative series in 3d gauge 
theories, considered as an effective theory yields a non- 



*In [||] it was shown that the convergence of the perturbative 
result for the free energy density can be improved using Pad'e 
approximants 



vanishing magnetic mass [j8|— |Tl]l . The 3d effective theory 
emerges from the high temperature limit of SU(N) gauge 
theory as follows: If the temperature is high enough the 
asymptotic freedom ensures the separation of different 
mass scales 2irT ^> gT ^> g 2 T and one can integrate 
out the heavy modes, with wave numbers ~ 2-kT and 
~ gT. This yields effective theories describing IR dynam- 
ics at scale gT (adjoint Higgs model) and g 2 T (3d pure 
gauge theory), respectively jL2j. The usefulness of the 
effective theory approach for unrevealing the source of 
breakdown of perturbation theory and solving the prob- 
lem of the perturbative IR catastrophe was illustrated 
in ||]. The effective theory approach was used to study 
non-perturbative correction to the Debye screening mass 
]l3| p^[ . The mass gap of pure 3d gauge theory can be 
related to the magnetic mass of hot SU (N) gauge theory 
through the dimensional reduction. The presence of the 
mass gap in 3d gauge theory needs some comments. In 
p3| it has been claimed that the magnetic mass cannot 
exist in a confining theory, since it would imply a perime- 
ter law for large Wilson loops. However, this argument 
is valid only at tree level and non-perturbative analy- 
sis of the 2+1 dimensional gauge theory Jl7],[l8| shows 
that both magnetic mass and non-zero string tension are 
present in the theory at the same time. 

Although progress has been made in understanding 
the high temperature dynamics in the case of QCD the 
usefullness of the effective theory approach is question- 
able, because the coupling constant is large for any realis- 
tic temperature and therefore the separation of different 
length scales is not obvious. 

In this contribution we will try to clarify whether 
screening masses defined as poles of the corresponding 
propagators can be determined in the 3d adjoint Higgs 
model, considered as an effective theory of hot SU(N) 
theory (e.g. QCD). We are going to review recent papers 
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of the authors |l^j2(| . The investigation of the screening 
masses defined as poles of the corresponding propagators 
is of great interest because they present essential input 
into the perturbative calculation of the thermodynamical 
quantities Jl]||. In it was shown that the contribu- 
tion of the magnetostatic sector to the free energy density 
can be understood in terms of massive quasiparticles with 
mass ~ g 2 T . The outcome of the resummed perturba- 
tive calculation for the free energy may depend on the 
correct choice of the screening mass [p2 23 . Therefore 
we will also discuss how far the values of the screening 
masses determined here influence the result of the per- 
turbative free energy calculation. We will proceed as fol- 
lows: in section 2 coupled gap equations for the screening 
masses are introduced and analyzed numerically. In sec- 
tion 3 the lattice formulation of the adjoint Higgs model 
is considered, numerical results for SU (2) Landau gauge 
propagators are discussed and compared with recent re- 
sults from 4d simulation of hot SU (2) theory. In section 
4 propagators are studied in the so-called A-gauges 
and numerical arguments for gauge independence of the 
pole mass within this class of gauges will be presented. 
In section 5 application of the pole masses to the resum- 
mation of the free energy is discussed. 



schemes are known to provide real values for the mag- 
netic mass fic|| . The corresponding expression for the 
on-shell self-energy reads 



IIx(fc = itot,tot) = CrriT, 
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Since we are interested in calculating the screening 
masses in the three-dimensional SU(N) adjoint Higgs 
model, Ht(/c,tot) should be supplemented by the corre- 
sponding contribution coming from Aq fields. This con- 
tribution is calculated from diagrams: 



and its analytic expression is 
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II. COUPLED GAP EQUATION 

As was discussed in the introduction a gauge in- 
variant definition of screening masses through the pole 
of the gluon propagator is possible and they can be 
determined self-consistently through the gap equations 
P,p|-pd]|. However, the determination of the electric and 
magnetic masses was attempted independently form each 
other. In view of the fact that g ~ 1 it seems natural to 
determine the screening masses from a coupled system of 
gap equations which does not rely on the separation of 
the electric and magnetic scales. We will assume, how- 
ever, the decoupling of non-static modes and the coupled 
gap equation will be derived from the effective theory, the 
3d adjoint Higgs model. The lagrangian of this theory in 
the Euclidian formulation is the following: 
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Gauge invariant approaches for the magnetic mass gen- 
eration in three-dimensional pure SU (N) gauge theory 
were suggested by Buchmuller and Philipsen (BP) M and 
by Alexanian and Nair (AN) [||. The approach of AN 
uses a hard thermal loop inspired effective action for the 
resummation in the magnetic sector. The approach of 
BP makes use of a gauged cr-model, and goes over to the 
SU(N) gauge theory in the limit of an infinitely heavy 
scalar field. At present only these two gauge invariant 



It is transverse and gauge parameter independent, it also 
does not depend on the specific resummation scheme ap- 
plied to the magnetostatic sector. It should be also no- 
ticed that it starts to contribute to the gap equation at 
0(g 5 ) level in the weak coupling regime, thus preserving 
the magnetic mass scale to be of order g 2 T. This is the 
reason why no " hierarchy" problem arises in this case in 
the weak coupling region. 

At 1-loop order the following diagrams contribute to 
the ^4()-self energy [] 




The self energy of Ao depends on the specific resum- 
mation scheme. For BP resummation it reads 
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In the AN scheme the corresponding expression reads 



' There is also a diagram arising from quartic self coupling of 
Aq, however, since the corresponding coupling is very small, 
we have not taken it into account. 
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FIG. 1. The temperature dependence of the scaled Debye 
mass for BP resummation scheme (solid) and for the AN re- 
summation scheme (dashed). The scaling factor is moo- 
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Although the value of IIoo is different for the two resum- 
mation scheme the analytic properties and on-shell values 
are the same. The coupled system of gap equations can 
be written as 
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The temperature dependence of mjj obtained from (g) 
is shown in Figure 1 for both schemes, where we have 
normalized the Debye mass by the leading order result, 
m do ■ The temperature dependence of the magnetic mass 
is shown in Figure 2, where we have normalized mx by 
the value of the magnetic mass obtained for pure three- 
dimensional SU (2) theory, in the BP (AN) gauge invari- 
ant calculations |M . As one can see the influence of Aq 
on the magnetic mass is between 1 and 3%. From Figures 
1 and 2 it is also seen that the temperature dependence 
of the screening masses is very similar to the temperature 
dependence of the respective leading order results. 
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FIG. 2. The temperature dependence of the scaled mag- 
netic mass for BP resummation scheme (solid) and for the 
AN resummation scheme (dashed). The scaling factors are 
m.y P and m^ N , respectively. 



III. LATTICE FORMULATION OF 3D ADJOINT 
HIGGS MODEL AND ITS PHYSICAL PHASE 

The lattice action for the 3d adjoint Higgs model used 
in the present paper has the form 

p L 

^^o(x)C/i(x)A (x + S)E/-t(x) + 



,(9) 



where Up is the plaquette, Ui are the usual link vari- 
ables and the adjoint Higgs field is parameterized by anti- 
hermitian matrices Aq = i J2 n a a A^ (a a are the usual 
Pauli matricies), as in ]T^ , |25| |T Furthermore /3 is the lat- 
tice gauge coupling, x parameterizes the quartic self cou- 
pling of the Higgs field and h denotes the bare squared 
mass of the Higgs particle. This model is known to have 
two phases the symmetric (confinement) and the bro- 
ken (Higgs) phase Q separated by the line of I s * order 
phase transition for small x. The strength of the transi- 
tion is decreasing as x increases and turns into a smooth 
crossover at x ~ 0.3 14 2%. 
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The high temperature phase of the SU(2) gauge the- 
ory corresponds to some surface in the parameter space 
(f3,x,h), the surface of 4d physics. This surface may lie 
in the symmetric phase or in the broken phase, i.e. the 
physical phase in principle might be either the symmet- 
ric or the broken phase. In the previous section it was 
tacitly assumed that the symmetric phase is the physi- 
cal one. This seems to be reasonable because m 2 DQ > 
and at tree level no symmetry breaking occurs. How- 



ever, it was shown that the 1-loop |28| and the 2-loop 
]l4[ effective potential of Aq has a non-trivial minimum, 
i.e. the broken phase might be the physical. The conclu- 
sion that the high-T physics is mapped onto the broken 
phase of the 3d model is self-contradictory has been ar- 
gued in a recent numerical lattice investigation, where 
the parameters appearing in ( |9|) w ere determined from a 
2-loop dimensional reduction [lj]. The dimensional re- 
duction performed in on the other hand led to the 
conclusion that the physical phase is the symmetric one. 
The expectation value of the adjoint Higgs field in the 
broken phase is of order ~ ^(g)' therefore for g << 1 
dimensional reduction is valid only if g?,AQ «T . What 
happens at intermediate couplings, however, remains an 
open question. 

We will discuss two possibility for fixing the param- 
eters appearing in (^|) by matching non-perturbatively 
some quantities which are equally well calculable both in 
the full 4d lattice theory and in the effective 3d lattice 
theory. We propose to match Landau gauge correlators 
of static link configurations calculated in the 4d theory 
to the results of the effective model. These are quantities 
calculated in gauges fixed independently. The assump- 
tion is that static configurations saturate also the corre- 
lators of the full theory. For such configurations the 4d 
Landau gauge is equivalent to its 3d version. 

In general this strategy requires the realisation of 
a matching in a three dimensional parameter space 
(/3, x, h), which is clearly a difficult task. We followed 
a more modest approach and fixed two of the three pa- 
rameters, namely (3 and x at values obtained from the 
perturbative dimensional reduction. The values of these 
parameters at 2-loop level arc jlj] 
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with a and T denoting the lattice spacing and the tem- 
perature, respectively. The coupling constant of the 4d 
theory g 2 (n) is defined through the 2-loop expression 
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For the comparison of the results of the 3d and 4d sim- 
ulations with it is necessary to fix the renormalization 
and the temperature scale. We choose the renormaliza- 
tion scale /i = 2t:T, which ensures that corrections to the 
leading order results for the parameters g\ and x of the 
effective theory are small. Furthermore we use the rela- 
tion T c = 1.06A-rjg from (7). Now the temperature scale 
is fixed completely and the physical temperature may be 
varied by varying the parameter x. 

Let us discuss the selection of the value of h for fixed (3 
(lattice spacing), i.e. the ways how to choose the line of 
4d physics. In our investigations three alternative choices 
for the line of 4d physics h(x) have been explored. A 
comparison with 4d simulations then allows the determi- 
nation of the physical line h(x). The first alternative for 
h(x) is the perturbative line of 4d physics, calculated in 
]l4| and lying in the broken phase, the other two are in 
the symmetric phase. The three alternatives are illus- 
trated on Figure 3 for (3 = 16. The actual procedure we 
used to choose the parameter h in the symmetric phase 
is the following. First, we determined the transition line 
htr{x). The transition line as function of x in the infinite 
volume limit was found in |14| ] in terms of the renormal- 
ized mass parameter y — rrc j 1 g\ (m is the continuum 
renormalized mass). It turns out to be independent of (3. 
Then using eq. (5.7) from @ one can calculate ht r {x). 
The use of the infinite volume result for the transition 
line seems to be justified because most of our simula- 
tions were done on a larger 32 2 x 64 lattice. The two sets 
of h{x) values, which appear on Figure 3, were chosen 
so that the renormalized mass parameter y (calculated 
using eq. (5.7) of ) always stays 10% and 25% away 
from the transition line. These values of h are of course 
ad hoc and one should use them only as trial values. 

Let us first discuss our calculations in the broken 
phase. The simulations were done for two sets of pa- 
rameters: (3 = 16, x = 0.03, h = -0.2181 and 
(3 = 8 , x = 0.09, h = —0.5159, here h was chosen along 
the perturbative line of 4d physics [JlJJ . The propagators 
obtained by us in the broken phase show a behaviour 
which is very different from that in the symmetric phase 
and that in the 4d case studied in Ref. The mag- 

netic mass extracted from the gauge field propagators is 
0.104(20)g| for the first set and 0.094(8)g| for the sec- 
ond set of the parameters. It is by a factor 4 to 5 smaller 
than the corresponding 4d result. Moreover, the prop- 
agator of the Aq field does not seem to show a simple 
exponential behaviour, what is in qualitative agreement 
with the findings of Ref. Q . These facts suggest that the 
broken phase can not correspond to the physical phase. 

We turn to the discussion of our results in the symmet- 
ric phase. In order to find the parameter range of interest 
for h we first have analyzed three different values of h at 
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FIG. 3. The bare masses normalized by the critical mass 
htr(x) used in our matching analysis (squares): the perturba- 
tive line (a) and two alternative sets / (line b) and (line c) . 
For a discussion of their choice see text. For x = 0.07 also a 
point deeper in the symmetric phase has been explored. The 
thick solid line is the transition line and its thickness indicates 
the uncertainty in its value. 

j3 = 16 and x — 0.07. The location of these values rela- 
tive to the transition line is shown in Figure 3. For the 
electric screening mass we find, in increasing values of h, 
m D /T = 1.72(10), 1.90(7) and 2.41(11). These results 
should be compared to the 4d data. From the fit given in 
Ref. 0, we find at T/T c — 12.57 for the electric screening 
mass mo/T = 1.85. This shows that our third value for 
h clearly is inconsistent with the 4d result. From a linear 
interpolation between the results at the three different 
values for h wc find the best matching value, i.e. a point 
on the line of 4d physics, h(x = 0.07) = -0.2496. 

The temperature dependence of the screening masses 
obtained in the symmetric phase for these two sets of pa- 
rameters, which stay close to the transition line, is shown 
in Figure 4. Also given there is the result of the 4d sim- 
ulations §, m 2 D /T 2 = Ag 2 (T), with A = 1.70(2) for the 
electric mass and m T /T = Cg 2 {T), with C = 0.456(6) 
for the magnetic mass. As can be seen both masses can 
be described consistently with a unique choice of the cou- 
pling h for temperatures larger than 10T C . Although even 
at T ~ 4T C we find reasonable agreement with the 4d fits, 
we note that the accurate choice of h becomes important 
and a simultaneous matching of the electric and magnetic 
masses seems to be difficult. For larger temperatures we 
find that the magnetic mass shows little dependence on h 
(in the narrow range we have analyzed) and the determi- 
nation of the correct choice of h thus is mainly controlled 
by the variation of the electric mass with h. 

Let us summarize our findings for the screening masses 
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FIG. 4. The screening masses in units of the temperature. 
Shown are the Debye mass tud for the first (filled circles) and 
the second (open circles) set of h, and the magnetic mass tut 
for the first (filled squares) and the second (open squares) set 
of h. The line (a) and line (b) represent fits to the temper- 
ature dependence of the Debye and the magnetic mass from 
4d simulations of [6]. The magnetic mass for the set II at 
temperature T ~ 90T C and ~ 9000T C was shifted in the tem- 
perature scale for better visualization. The open triangle is 
the value of the Debye mass for x — 0.07 and h — —0.2179. 

in the symmetric phase. For T > 10T C the screening 
masses can be described very well in the effective theory. 
Suitable values of h can be found using the interpolation 
procedure outlined above for x = 0.07. This procedure 
can also be followed for x = 0.05 and 0.03. There the 
4d data are well matched for values corresponding to set 
I (see Figure 4), therefore the following h values can be 
considered as ones corresponding to 4d physics, h{x = 
0.07) = -0.2496, h(x = 0.05) = -0.2365 and h(x = 
0.03) = —0.2085. An interpolation between these values 
gives the approximate line of 4d physics ti4,d(x). 

IV. GAUGE INDEPENDENCE OF THE 
PROPAGATOR POLE MASS 

In confining theories the propagator pole does not cor- 
respond to an asymptotic state (unlike e.g. in QED) 
therefore there is no physical reason for gauge indepen- 
dence of the pole mass. However, the propagator pole 
was proven in the framework of perturbation theiry to 
be gauge independent in the high temperature deconfincd 
phase of SU(N) theory |§. 

We have investigated the gauge dependence of the ef- 
fective masses numerically using the so-called A-gauges 
pi| , defined by the gauge fixing condition 
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Xd 3 A 3 + d 2 A 2 + d x A x = 0. 



(14) 



Case A = 1 corresponds to the Landau gauge. The possi- 
ble gauge dependence of the pole mass has been checked 
by investigating in addition to the Landau gauge prop- 
agators also the propagators for A = 0.5 and 2.0 for 
/3 = 16, x = 0.03 and h = -0.2085 on a 32 2 x 96 
lattice. The results of these measurements together with 
the results obtained from Landau gauge propagators (on 
32 2 x 64 lattice) are summarized in Table 1. 

/3 = 16 



A 


0.5 


1.0 


2.0 


m D /g'i 


1.22(5) 


1.32(5) 


1.20(10) 


m T /g'i 


0.48(5) 


0.46(8) 


0.42(5) 



TAB 1: Numerical investigation of the dependence of the 
pole mass on the gauge parameter A. 



V. FREE ENERGY RESUMMATION 

In this final section we will use the screening masses 
calculated before for the evaluation of the free energy 
density. The partition function can be calculated in the 
effective theory in the following way: 
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where Z non ^ stat is the contribution of the non-static 
modes to the partition function which was calculated in 
H to 3-loop level, and L e ff is the lagrangian of the ef- 
fective theory given by (1). Performing the integration 
over static electric fields (Aq) one obtains the contribu- 
tion of the electric scale (gT) to the partition function 
Z e i. This step can be performed perturbatively because 
Aq has a non-zero thermal mass. The resulting contribu- 
tion was calculated in || yielding odd powers to the weak 
coupling expansion of InZ. Integration over static mag- 
netic fields yields the contribution of length scale g 2 T 
to the partition function and was calculated using nu- 
merical lattice simulation in |2l| . In this way the weak 
coupling expansion can be obtained to 0(g e ). However, 
there are at least two things one may worry about in the 
outlined perturbative program : i)In the calculation of 
Z e i the tree-level (from the point of view of the effective 
theory) mass mno was used. As we have seen both lattice 
simulations and coupled gap equations yield a mass for 
the Ao field, which is considerably larger than the tree- 
level result rriDo- ii) For g ~ 1 the separation of electric 
(gT) and magnetic scales (g 2 T) is not obvious and it is 
interesting to investigate the influence of the latter on 
the A integration. This investigation is also motivated 
by the fact that the electric mass is rather sensitive to 
the magnetic scale. To investigate the effect of using the 
'exact' masses we will reorganize the perturbation theory 



and perform a loop expansion, instead of the expansion 
in powers of g. This can be achieved by rewriting the 
Lagrangian as 
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are the 2- and the 3-loop counterterms to be treated as 
perturbations. The values of mo and mj will be taken 
from the coupled gap equation using the BP scheme. At 
2-loop level the following diagrams will contribute to the 
free energy, — h\Z e i/V, 






(The solid line denotes the Aq propagator, the dot de- 
notes the Aq — Ai vertex and the cross denotes the mass 
counterterm) . We will also study the 3-loop contribution 
in the limiting case m<r = 0. In this case the following 
3-loop diagrams should be taken into account in addition 
to those calculated in 11: 




The contribution of these diagrams to f e i 
are: 



r 



(N 2 -l)Ngl (0 2 l m 2 



2mjjmT + (4TO 2 } — 
f b = - 

f c = - 



ln 



A 



'" - - - kN i J)f g3 m D m T 

c _ (A M -l)(mp -roj)ro_ D 

fi = -\ A (N*-l)^ 



■hxZ el /V 

(20) 
(21) 
(22) 
(23) 



6 




500 



T/T,. 



FIG. 5. The resummed loop expansion of the SU(3) free 
energy: a) the 1-loop, b) 2-loop and c) the 3-loop level free 
energy for hit = 0. 
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FIG. 6. The weak coupling expansion of the SU(3) free 
energy: a, b, c and d are the 0(g 2 ), 0(g 3 ), 0(g 4 ) and 0(g 5 ) 
level results 



which shows the known alternating behaviour in different 
orders. Compared to this in the resummed loop expan- 
sion the contribution of the static modes is larger, but the 
free energy decreases systematically, contrary to the al- 
ternating behaviour of the weak coupling expansion. The 
3-loop level resummed free energy if compared to 0(g 5 ) 
order result differ from it by the amount of a few percent. 
Finally the effect of massive magnetic modes was studied 
at 2-loop level and it turns out that the variation of the 
free energy is about 1% relative to the tot = case. 



VI. SUMMARY 

The free energy density has been shown in this paper 
to be rather sensitive to the actual values of the screen- 
ing masses. This is probably also true for other physical 
quantities. This gives the motivation for the careful an- 
alytical and numerical investigation of these quantities, 
presented in our paper. 

Both numerical Monte-Carlo simulations and investi- 
gations of the coupled gap equations show that the mag- 
netic mass of the 3d adjoint Higgs model is very close to 
the magnetic mass of the pure 3d gauge theory, however, 
the numerical values of the magnetic mass obtained in 
these aproaches are different. The 1-loop gap equation 
approach gives the magnetic mass roughly equal to 0.2 853 
for BP scheme and 0. 38173 f° r AN scheme, while the nu- 
merical lattice simulations give 0.46g|. The temperature 
dependence of the Debye mass was found very close to the 
leading order result both in the gap equation approach 
and in the numerical lattice simulations, but the numeri- 
cal values are again different and equal to (1.2 — 1.3)muo 
for the coupled gap equations ( for T > 100T C ) depend- 
ing on resummation scheme (see Figure 1) and 1.6touq 
for lattice simulations. 

We leave the interpretation of these systematic devia- 
tions to a future publication. 
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First let us investigate the resummed loop expansion for 
the free energy in the case tot = 0. The numerical re- 
sults of the loop expansion are shown in Figure 5. For 
comparison we have also plotted in Figure 6 the weak 
coupling expansion for the free energy calculated in 
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